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ABSTRACT 

We report an SMA interferometric identification of a bright submillimeter source, GOODS 850-5. This 
source is one of the brightest 850 /im sources in the GOODS-N but is extremely faint at all other wavelengths. 
It is not detected in the GOODS HST ACS images and only shows a weak 2 <t signal at L4 GHz. It is detected 
in the Spitzer IRAC bands and the MIPS 24 fim band, however, with very low fluxes. We present evidence in 
the radio, submillimeter, mid-IR, near-IR, and optical that suggest GOODS 850-5 may be a z > 4 galaxy. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies : star- 
burst — IR: galaxies — submillimeter 



1. EVITRODUCTION 

The Rayleigh- Jeans portion of the dust spectral energy dis- 
tribution (SED) produces a strong negative /T-correction and 
makes the observed submillime ter flux of a dusty gala xy al- 
most invariant atz> Itoz^lO jBlain & Longairll993b . This 
makes the submillimeter wavelength a potentially powerful 
probe to the high-redshift universe. However, to date, all the 
identified submillimeter galaxies are at redshifts lower than 
4, likely because of the limited resolution of current submil- 
limeter instruments and the limited sensitivity of current radio 
instruments. 

The Submillimeter Common-User Bolometer Array 
(SCUBA) on the single-dish James Clerk Maxwell Tele- 
scope resolved 20%-30% of the submillimeter extragalactic 
background light into point sources brighter than ^ 2 mJy 
at 850 i-im (Barger, Cowie, & Sanders 1999; Eale s et aTl 
Il999l) . Because of the low resolution (- 15") of SCUBA, 
identifications of the submillimeter sources have to assume 
the ra dio-FIR correlation in local galaxies (see, e.g.. lCondonl 
Il992h and rely on radio interferometry to pinpoint the 
location of the submillimeter emission. Optical spectroscopy 
of radio identified submillimeter sources shows that they 
are ultraluminous (> 10'^ Lq, corresponding to star for- 
mation rates of lO^-lO-' M© yr"') sources at z ~ 2-3 and 
that they d ominate the total star formation at this redshift 
d Chapman et al. 2005). However, the positive /T-correction of 
the radio synchrotron emission makes the radio wavelength 
insensitive to high-redshift galaxies, and radio observations 
can only identify 60%-70% of the blank-field submillimeter 
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sourc es (e.g., iBarger. Cowie. & Richardsl [20001 llvison et aP 
I2OO2I). The radio unidentified submillimeter sources are 
commonly thought to be at redshifts higher than the radio 
detection limit (typically z ^ 3-4) but there is no direct 
evidence for such a high-redshift radio-faint submillimeter 
population. 

With recent development in submillimeter interferometry, 
it is now possible to locate submillimeter sources directly 
without relying on radio observations. We have begun a 
program to target radio-faint submillimeter sources with the 
Submillimeter Array (SMA^) to determine whether there is 
a high-redshift (z > 4) tail in the redshift distribution of 
the submillimeter sources. Here we report our first iden- 
tification in this program, GOODS 850-5. GOODS 850- 
5 was detected by our SCUBA jiggle-map survey in the 
Gr eat Observatorie s Origins Deep Surveys-North (GOODS- 
N,'Giavalisco et al."2004a) with an 850 ^m flux of 12.9 ±2.1 
mJy (Wang, Cowie, & Barger 2004). It was also detected in 
the combined jiggle and scan map of GOODS-N (GN 10, 
see IP ope et al.l l2006l and references therein). It is the sec- 
ond brightest submilUmeter source in our jiggle-map catalog 
of the GOODS-N and has an IR luminosity of ~ 2 x 10'^ 
Lq. It does not have a 5(t radio counter part in the deep Very 
Large Array (VLA) 1. 4 GHz catalogs of Richardsl (l2000l) and 
iBiggs & IvisonI (l2006i) . We obtained an unambiguous iden- 
tification of GOODS 850-5 with the SMA and found that 
the counterpart to this source is remarkably faint at all other 
wavelengths. All data in the optical, near-IR, mid-IR, submil- 
limeter, and radio point to a source at z > 4. This is impor- 
tant evidence for the existence of high-redshift submillimeter 
sources. In this letter, we report the SMA observation and 
the multiwavelength photometry (§ |2]i and the redshift con- 
straints derived from the photometric data (§ |3ll. We discuss 
the implication for the star formation history in §|4]and sum- 
marize in §|5] We adopt Hq = 71 km s"' Mpc"', = 0.73, 
and = 0.27. 

2. OBSERVATIONS AND MEASUREMENTS 

A full track of SMA observations were obtained on 2007 
January 23 using eight antennas in the compact configura- 
tion. The receivers were tuned to a frequency of 345.0 GHz 
centered in the upper sideband, yielding 335.0 GHz in the 
lower sideband. Each sideband covers a 2 GHz bandwidth 
and contains 768 spectral channels. The tracking center 
is at the SCUBA position of GOODS 850-5 (a[J2000.0] = 
12''36'"33\45, (5[J2000.0] = 62°14'09"43). Callisto and 
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Fig. 1 .— Multi-wavelength images of GOODS 850-5. Each panel has a size of 24". North is up. All four of the ACS (F435W, b; F606W, v; F775W, i; and 
F850LP, z) bands and all four of the IRAC channels (1, 2, 3, and 4) are included in the color pictures. The color codes are labeled in the images. Grayscale 
images have inverse scales. The SMA position is labeled with 2" diameter circles. The contours in the SMA image have levels of -2, 2, 4, and 6 cr where 1 cr is 
1.4 mjy beam"' . 

the bright quasar 3C84 were observed to provide flux and 
bandpass caHbrations, respectively. Quasars 1419-1-543 and 
1048-1-717, respectively 15?5 and 14° away from the target, 
were observed after every 15 minutes of on-target integration 
for time-dependent complex gain calibrations. The 225 GHz 
opacity was an excellent 0.03-0.06, and the dual-sideband 
system temperature averaged over the entire track was 188 
K. 

The caUbration and data inspection were performed with 
the Caltech package MIR, modified for the SMA. Continuum 
data were generated by averaging the spectral channels after 
the passband calibration. Both gain calibrators were used to 
derive gain curves. However, we carried out a separate reduc- 
tion in which the primary calibrator 1048-1-717 was used to 
calibrate 1419H-543. The reduced image of 1419H-543 has a 
point source 0."24 away from the phase center. Since the two 
calibrators are separated by a large 28° on the sky, this offset 
should be only considered as an upper limit to the astrometry 
error in the reduced image of GOODS 850-5. The calibrated 
visibility data were Fourier transformed and deconvolved in 
the package MIRIAD to form images. In the transformation, 
we weighted each visibility point inversely proportional to the 
system temperature. Combined with a natural weighting, this 
provides an optimum S/N and a synthesized beam FWHM of 
2 "4 X 2. "2 at a position angle of 24°. In the dirty image of 
GOODS 850-5, a point source is apparent and we restricted 
the CLEAN deconvolution to a ^ 10" region centered at the 
point source. There is no sign of other sources both in the 
dirty and cleaned images. The noise level in the cleaned im- 
age is 1 .4 mJy beam"' . 

A point source of 12.0 ± 1.4 mJy is clearly detected very 
close to the phase center. A Gaussian fit to the image yields a 
source profile identical to the synthesized beam, indicating an 
unresolved source with good calibrations. The fit provides a 
source position of a(J2000.0) = 12''36™33^45, (5(J2000.0) = 
62°14'08"65, which agrees extremely well with the SCUBA 
position. The positional error in the fit is 0"2, consistent with 
the angular resolution for this S/N and also comparable to the 
aforementioned upper limit to the astrometry error derived 
from the calibrator 1419-1-543. We also performed a point- 
source fit directly to the visibility data, and the result is fully 
consistent with the fit in the image plane. 

We present the GO ODS HST Advanced Ca mera for Sur- 
veys (ACS) images dGiavalisco et al.l l2004ab . our CFHT 
WIRCam Ks band image (36 hours of integration, Keenan et 
al. 2007, in preparation), the GOODS Spitzer Legacy Program 
(Dickinson et al. 2007, in preparation) Infrared Array Cam- 
era (IRAC) and Multiband Imaging Photometer for Spitzer 
(MIPS) images, the VLA image (Biggs & Ivison 2006), and 
our SMA image of GOODS 850-5 in Fig.[T] In all of the four 
ACS images (F435W, F606W, F775W, and F850LP) and the 
WIRCam image, there are no galaxies within 0."5 from the 



SMA position. A faint galaxy is clearly detected by Spitzer 
in all four IRAC bands and in the 24 ^m MIPS band right 
at the SMA position. Because of the high accuracy (0."2) of 
the SMA astrometry, this Spitzer source is a robust identi- 
fication of GOODS 850-5. We list our broadband fluxes in 
Table [U The 850 /im flux is a noise weighted mean from the 
SCUBA jiggle map (Wang et al. 2004) and the SMA measure- 
ment. Typical errors in flux calibrations at this wavelength are 
> 10% and could be larger than the noise. However, the fact 
that the SMA flux and the SCUBA flux agree with each other 
extremely well suggests good flux calibrations in both exper- 
iments. The 1.4 GHz flu x is measured at the SMA position 
from the VLA image of Biggs & Ivison (2006), assuming a 
point source and with a primary beam correction. The ACS 
and WIRCam fluxes are also measured at the SMA position 
with respectively 0."8 and 1"5 diameter apertures. Measur- 
ing fluxes in the Spitzer bands is more difficult. In the 3.6 
and 4.5 /im images, there is a nearby blue galaxy blended 
with GOODS 850-5. There is another galaxy blended with 
GOODS 850-5 in the MIPS image. To deconvolve the fluxes, 
we performed PSF fitting to the blended sources. We used 
the IRAC and MIPS in-flight PSFs and applied simultaneous 
two-component fits to the two sources, and scaled the mea- 
sured fluxes to total fluxes using the PSFs. To verify the fluxes 
measured in this way, we first compared them with simple 
aperture-corrected aperture fluxes. We found that at 5.8 and 
8.0 urn where GOODS 850-5 is isolated, the PSF fitting pro- 
vides fluxes consistent with the aperture fluxes. We also used 
the SExtractor package (Bertin & Arnouts 1996) to measure 
fluxes. SExtractor could deblend the two sources at 4.5 /im 
but not at 3.6 /im even with the most aggressive deblending. 
The aperture corrected SExtractor fluxes at 4.5 /im and longer 
are all consistent with the PSF fitting fluxes. We thus conclude 
that our PSF fitting provides reliable estimates of the fluxes of 
GOODS 850-5. The Spitzer fluxes Hsted in Table [T]ai-e all 
measured with the PSF fitting method. 

We note that o ur SMA identificat ion confirms the counter- 
part suggested bv lPope et al.l (l2006h . Our carefully measured 
Spitzer fluxes are also consistent with the SExtractor pipeline 
fluxes in Pope et al. The 1.4 GHz flux in Pope et al. is 76% 
larger than ours but the difference is within the errors. 

3. REDSHIFT CONSTRAINTS 

GOODS 850-5, one of the brightest 850 /im sources in the 
GOODS-N, is incredibly faint at all other wavelengths. Be- 
cause of the optical and near-IR faintness, it is difficult to 
obtain its redshift through spectroscopy. We obtained red- 
shift constraints for GOODS 850-5 using the multiwavelength 
broadband fluxes. There are three independent pieces of evi- 
dence that suggest GOODS 850-5 is a z > 4 galaxy. 

By assuming that the well known radio-FIR correlation 
for local star forming galaxies also applies to submillime- 
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TABLE 1 
Photometry of GOODS 850-5 



'5f435W 




Sfiisw 


5f850LP 




^3.6^111 




'^5.8^111 






S&SQ^im 5l.4GHz 


-0.0 13 ±0.004 


-0.004±0.003 


0.001 ±0.006 


-0.009±0.009 


0.055 ±0.052 


1.14±0.14 


1.64±0.13 


2.33±0.24 


5.37±0.37 


46.3±9.2 


12300±1200 18.7±8.0 



Note. — All fluxes are in fiJy, which can be converted to AB magnitude with Afi = 23.9- 2.5 log(/iJy). 



ter sources, we can estima te the redshift usin g the radio and 
submillimeter fluxes (e.g., ICariUi & Yunlll99 9: Barger et al. 
2000). Barger et al. assumed pure power-law SEDs in both 
the radio and the submillimeter and obtained the formula 
z = 0.98(5850,.m/5i.4GHz)" -''- 1, based on Arp 220. This im- 
pUes z ~ 4.3 for GOODS 850-5. This is a lower limit, because 
the peak of the dust SED and the free-free emission start to 
move in the observed bands at z > 4. Both effects underes- 
timate the redshift. By ad opting the complete Arp 220 SED 
model in lSilva et al.1 (fT998), we find z - 5 .6 for G OODS 850- 
5 (z ~ 4.3 if we adopt the radio flux of 33. 1 /iJy in lPope et al.l 
l2006h . This redshift estimate is nevertheless very crude. Be- 
cause of the uncertainties in the submillimeter and radio SEDs 
and the radio-FIR correlation at high redshift, s uch redshift 
estimates typically have an error of ±0.5 (e.g., Ilvison et"an 
l2002h . 
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Fig. 2.— Sgsoa m/S24um Aux ratio vs. redshift. Curves are filter-convolved 
ISilva et al] 419981) model SEDs of Arp 220 (ultraluminous starburst, solid 
curve), M 100 (normal spiral, dash-dotted cun>e), and M 82 (low luminos- 
ity starburst, dashed ciin'e), and Mrk 231 (ultraluminous dusty AGN, dotted 
cun'e) derived from the photometry in NED** . Shaded region is the range of 
the models i mCharv & Elbazi t2001.) . Solid squares are GOODS-N S CUBA 
sources identified at 24 fim or 1.4 GHz (Wang. Cowie. & BargeJ200^ . Hor- 
izontal lines are GOODS 850-5 and the 1 a error. 

The unusually large Sfis0fj,m/S24f_im ratio also suggests a high 
redshift of z > 3. As can be seen in Fig. |2l the Ss.5o^m/S24fj.m 
ratio of GOODS 850-5 is an order of magnitude higher than 
the mean value for other identified submillimeter sources in 
the GOODS-N. This can only be explained by a z > 3 galaxy 
with any dust SEDs allowed by current models. 

Finally, we used the photometric redshift technique to test if 
the HST, WIRCam, and Spitzer SED of GOODS 850-5 can be 
fitted with a z > 4 galaxy. Because there is not a unique cor- 
relation between the dust SED and stellar SED of galaxies, 
we only included the ACS, WIRCam, and IRAC passbands 
in the photometric redshift study to account for just the stel- 
lar emission. Although GOODS 850-5 is not detected in the 
ACS images, the upper limits to the ACS fluxes suggest a very 
red SED at < 3 fim (Fig. [3]) and thus still provide useful con- 
straints. In the photometric redshift fitting, we assigned zeros 
to the ACS fluxes and used the 1 cr eiTors as the flux uncer- 
tainties. We used the measured flux, even though it is only 
1 cr. 

We used the popular Hyperz package 
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Fig. 3. — Photometric redshift result for GOODS 850-5. The upper panel 
shows the observed IRAC, WIRCam, and ACS SEDs (squares), and the best- 
fit SEDs (solid curves) in the empirical set (red) and in the population syn- 
thesis set (blue). The lower panel shows the (dashed curves), the corre- 
sponding probability (solid curves), and the extinction {dotted curx'es) of the 
best-fit SEDs vs. redshift for the two SED sets (same color coding as in the 
upper panel). 



(Bolzo neUa. JVliraUes. & Pellol 12000^. We limit ed flie 
redshift range to z = 0-12 and adopted the Calzett Tet al.l 
(12000) extinction of Ay = 0-5. We used two sets of SEDs. 
The first set consists of empirically observed SEDs, including 
the Im, Scd, Sbc, and E type galaxies in C oleman et al. 
(1980) and the starburst SEDs SB2 and SB3 in lKinnev etall 
( 1996). The primary solution provided by Hyperz is z = 6.99, 
Ay = 0.5, and a Sbc type galaxy. The corresponding proba- 
bility is 65.1%, and the 99% confidence range is z = 5.6-8.0. 
The second SED set is the stellar population synthesis models 
include d in Hyperz, which is based on the earlier version of 
iBruzual & Chailot ( 1993). The primary solution is z = 6.42, 
Ay = 1.5, and a burst type SED with an age of 0.5 Gyr. The 
coiTesponding probability is 96.4%, and the 99% confidence 
range is z = 5.2-8.0. In both SED sets, the photometric 
redshift confidence range is consistent with the range allowed 
by the radio, submillimeter, and mid-IR flux ratios discussed 
above. In Fig.[3]we show the fitted and observed SEDs, and 
the best-fit probability, and extinction as functions of 
redshift. The SED can only be fitted with Ay > 2 at z < 3 and 
the probabilities at z < 3 of the best-fit models are zero even 
when Ay = 10 is allowed in the fitting. Therefore, a highly 
reddened galaxy at z < 3 is ruled out for GOODS 850-5. In 
the photometric redshift fitting, an intermediate redshift of 
3-4.5 is not entirely ruled out, although the probability is low 
(< 10%). 

The primary uncertainty in the above analyses is the pos- 
sibility of the existence of an active galactic nucleus (AGN), 
whose SED may be different then the ones used for the radio, 
submillimeter, IR, and optical photometric analyses. We ar- 
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gue that whether or not there is an AGN is not Hkely to change 
the high redshift of GOODS 850-5. First, GOODS 850-5 is 
not d etected in the 2 Ms Chandra images ([Alexander et alj 
l2003h in either the hard or soft X-ray bands. This rules out 
a Compton-thin AGN at z < 3. (The 2 Ms Chandra obser- 
vations can detect a lO**^ and 10"^^^ erg s"' soft and hard X- 
ray source, respectively, at z 3.) The best example of a 
local Compton-thick AGN in the far-IR luminosity class of 
GOODS 850-5 is Mrk 231. Although the IRAC power-law 
spectral slope of GOODS 850-5 is similar to that of Mrk 231, 
its S'i^Q^m/StAfj.m flux ratio is more than an order of magnitude 
higher than Mrk 231 atz < 3 (Fig.EJ, i.e., GOODS 850-5 does 
not show much of the warm dust emission that is characteris- 
tic for a heavily obscured AGN. Furthermore, the rest-frame 
1 .4 GHz flux of Mrk 23 1 also follows the radio-FIR corre- 
lation for normal galaxies with a q value of 2.24 (under the 
definition of q in ICondon|[l992l) but it has a flat radio spec- 
trum (a smaller A'-correction in the radio). Therefore, even if 
GOODS 850-5 is an AGN that is similar to Mrk 231, its large 
5'850/jm/5'i.4GHz flux ratio will suggest a redshift even higher 
than that estimated from the Arp 220 template. 

4. IMPLICATIONS FOR THE STAR FORMATION HISTORY 

There is growing evidence for a signific ation submillimeter 
population at z > 4. iDunlop et al.l (|2004|) presented evidence 
for an optically faint z ~ 4.1 galaxy for H DF 850.1. More 
recent SMA results of I Younger et al.l(l2007h suggest a number 
of radio and optically faint submillimeter sources to be at z > 
4. GOODS 850-5 is similar to these galaxies. On the other 
hand, GOODS 850-5 has the highest 585()^m/'Si .4GHz flux ratio, 
and also has higher quality data at many other wavelengths. 
Unlike HDF 850. 1, which is behind a bright galaxy, GOODS 
850-5 can be studied in details especially in the optical and IR 
wavelengths. 

The discovery of GOODS 850-5 and its like extends the 
redshift distribution of bright submillimeter sources to beyond 
4.0, and perhaps even to z > 5. The 12.3 mJy 850 jim flux 
of GOODS 850-5 implies an IR lum inosity of ^ 2.3 x 10'^ 
Lq (Lir = 1.9 X lO'^^ssoMini-o/mJy, iBlainet al.ll2002l) . and 
therefore a very high star formation rate of ^ 4000 Mq yr"' 
(M= 1.7 X 1 0" ' °L,p /Lg ■ iKennicut jl 1 998h . This is an extraor- 
dinary amount of star formation at high redshift but it can 
only be identified in the submillimeter Even if such a source 
could be picked up in optical and near-IR broadband dropout 
surveys, its star formation would still be missed due to the ex- 
treme faintness in the rest-frame UV. How much can hidden 
sources like this affect the current picture of the cosmic star 
formation history at high redshift? 

Bright submillimeter sources (> 2 mJy at 850 /Ltm) con- 



tribute 10 Jy deg"^ to the total background, and 60% to 
70% of them have radio counterparts. We assume that all 
of the ~ 4 Jy deg"^ radio-unidentified bright submillimeter 
sources are at z = 4-7. With the above conversions between 
850 /im flux, IR luminosity, and star formation rate, this im- 
plies a mean comoving IR luminosity density of ^ 2.2 x 10^ 
Lq Mpc"^, and therefore a star formation rate density (SFRD) 
of ^ 0.04 Mq yr"' Mpc"^. This is consistent with the es- 
timate in Barger et al. (2000) based on a smaller sample of 
radio unidentified submillimeter sources, after taking into ac- 
count the differences in cosmology. Our value of 0.04 M© 
yr~' Mpc"^ averaged over z = 4-7 is ~ 2 x higher than the typ- 
ical value of 0.01-0.02 Mq yr"' Mpc"^ af z ~ 4 determined 
from optical surveys without extinction correction (e.g., 
'Giava lisco et al.ll2004 b; Bouwens. Broadhurst. & lUingworthl 
[2003i) . In general, Lyman-break galaxies are not bright 
subm i llimeter sources dPeacock et al] 120001 : IChapman et alj 
120001: IWebb et al.ll2003h . suggesting that most (if not all) of 
the above SFRD estimated from the radio-unidentified bright 
submillimeter sources is missed by optical surveys. Although 
our SFRD of 0.04 Mq yr"' Mpc"^ is likely an upper limit due 
to the assumptions, it still shows the significance of dust hid- 
den star formation at high redshift. 

5. SUMMARY 

We obtained accurate astrometry for the bright submillime- 
ter source GOODS 850-5 with the SMA interferometer at 850 
/im. The counterpart is extremely faint in the optical, near- 
IR, mid-IR, and radio. Both the radio and 24 /im faintness 
of GOODS 850-5 suggests a high redshift of z > 4. The very 
red SED between 3.6 and 8.0 //m and the nondetection in the 
optical and Kg bands can only be fitted by stellar continuum at 
z > 4. This discovery provides important evidence that some 
of the radio undetected submillimeter sources are at high red- 
shift and extends the redshift distribution of bright submil- 
limeter sources to z > 4. It suggests that a great fraction of 
star formation at high redshift is hidden from optical surveys. 
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